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On  a Characterization  of  Multivariate  Distributions 
with  Applications  in  Reliability  and  Epidemiology 

by 

Naftali  A.  Langberg 
ABSTRACT 


Let  Tj,  . ..,  T be  positive  random  variables  with  finite  means.  Further 

let  I be  the  collection  of  all  subsets  of  {1,  . ..,  n},  and  let  £ be  a function 

from  the  nth  Euclidian  space  to  J,  that  equals  to  J,  (J  e I)  at  (a^,  . ..,  an) 

iff  min  a.T.  < min  a.T. . We  prove  that  min  a.T.  and  £(a.,  ....  a ) are  inde- 
ieJ  1 1 i^J  1 1 isisn  11  1 n 

pendent  random  variables  for  every  n real  numbers  aj , . . . , a^  iff  for  every  n 

positive  real  numbers  b^,  ....  bn  and  r = 1,  n the  random  variables 

min  a.T. /E(  min  a.T.)  and  T /ET  are  identically  distributed.  Further  we  provide 
lsisn  1 1 isisn  1 1 r r 

an  explicit  formula  for  the  distribution  of  £ (a.,  ...,  a ).  Multivariate 

i n 

distributions  that  possess  the  independence  property  are  presented.  Their  use 
in  Reliability  growth  or  decay  models  as  well  as  in  Mathematical  Epidemiology 
are  discussed. 


Key  Words:  Multivariate  distribution,  minima,  Weibull,  distribution,  independence 


1.  Introduction  and  Summary. 

Introduction.  Let  , ....  be  positive  random  variables,  and  let 
al’  '***  an  Positive  real  numbers.  If  T^,  . . . , are  the  initial  life 
lengths  of  n components  in  a series  system,  then  a^Tj,  ...,  can  be 

regarded  as  the  life  lengths  of  those  components  at  some  phase  of  a relia- 
bility growth  or  decay  process.  Let  an  individual  be  exposed  to  n contagious 
diseases  in  an  environment  consisting  of  infectives  and  susceptibles.  Then 
a.T, , ...,  a T may  describe  the  times  until  that  individual  becomes  an 
infective  from  disease  1 through  n respectively.  If  we  observe  a series 
system,  or  an  individual  who  is  exposed  to  n diseases,  only  two  quantities 
are  identifiable:  (i)  Time  until  occurence  (failure,  or  infection)  (ii) 

Cause  of  occurence  (failure  due  to  some  components,  infection  by  some  of 
the  diseases) . The  stochastic  representation  and  analysis  of  the  described 
models  simplifies  if  for  every  n positive  real  numbers  the  following  two 
properties  hold. 

Time  to  occurence  and  ca”se  of  occurence  are  independent  random  variables.  (1.1) 

" 

Time  to  occurence  and  cause  of  occurence  have  "identifiable"  distributions. (1.2) 

Recently  Langberg,  Lanzdorf  and  Proschan  (1978)  used  multivariate  distri- 
butions that  satisfy  (1.1)  and  (1.2)  to  describe  and  analyze  a variety  of 
reliability  growth  and  decay  models.  In  the  cited  reference  the  authors  con- 
sidered multivariate  distributions  with  independent  and  dependent  components, 
as  well  as  distributions  that  may  oi  may  not  be  absolute  continuous.  A well 
known  family  of  n-dimensior.al  random  vectors  with  independent  components  that 
satisfy  (1.1)  and  (1.2)  is  the  exponential  one.  More  specifically  let 

T, , ...»  T be  independent  exponential  random  variables  with  means  y. , ...,  y 
in  in 

respectively,  and  let  a^,  ...,  an  be  positive  real  numbers.  Then  rhe  following 
three  statements  hold. 

J 


L 


2 


41  111 

Time  to  occurance  is  exponential  with  mean  equal  to  ( £ a7  u-  )~  . 

j*l  J J 

T* 

-1  -1, 


-1  -K-l 


Tlie  probability  of  occurence  due  to  cause  i equals  to  (a  7 A p 7 )(  £ aTyTj 

11  j=l  3 3 

for  i * 1,  . ..,  n. 


Time  to  occurence  and  cause  of  occurence  are  independent  random  variables. 
Billard,  Lacayo  and  Langberg  (1978),  Lacayo,  Langberg  (1978)  and  Langberg 
(1978),  utilized  those  properties  of  independent  exponential  random  variables 
to  describe  and  analyze  n-dimensional  simple  epidemics. 

Summary : Lemmas  instrumental  to  the  proofs  of  our  main  results  are 
presented  in  Section  2.  In  Section  3 we  show  that  condition  (1.1)  is  equivalent 
to  the  following  statement. 

For  every  n positive  real  numbers  a.,  ....  a , (E  min  a.T.)-1  min  a.T.  and 

1 " lsisn  1 1 ma,  1 * 

(ETj)  *Tj  are  identically  distributed  random  variables  for  j = 1,  ...,  n. 

Clearly  (1.6)  provides  an  explicit  formula  for  the  distribution  of  the  random 
time  until  occurence.  In  addition  we  derive  in  Section  3 from  (1.6)  an  explicit 
form  for  the  distribution  of  the  cause  to  occurence.  Some  multivariate  distri- 
butions which  satisfy  (1.6)  are  presented  in  the  last  Section. 


(1.3) 


(1.4) 

(1.5) 


(1.6) 
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2.  Preliminaries. 

Let  denote  by  (T^,  T2)  a positive  random  vector  with  means  equal  to  p^,  p2 
respectively,  by  T(x)  the  minimum  of  x Tjt  and  T2>  and  by  F(x)  the  expected 
value  of  T(x).  Further  let  F'(x+)  and  F'(x-)  be  the  right  and  left  side 
derivatives  of  F at  the  point  x.  Finally  let  G(x)  be  equal  to  p^ICx  2:  O)F(x), 
Gn(x)  be  the  convolution  of  G and  a uniform  distribution  on  C—  07,  and  let 
gn(t)  be  the  density  function  of  G^t),  n = 1,  2,  ...  . For  reference  purposes 
we  summarize  in  Lemma  2.1  without  proofs  some  straight  forward  results. 


Lemma  2,1.  The  following  six  statements  hold. 

F is  a concave  r.ondecreasing  function. 

Lim  F (x)  = 0,  and  Lim  F(x)  = p 
x-K)+  x-x» 

F'(x+)  = ETj I (T2  > xTp  and  F:(x-)  = ET2I(T2  a xTj). 

Lim  Sup |G  (x)  - G(x)|  = 0. 
n-*»  x 

En  is  nondecreasing  in  n,  and  Lim  g (x)  = G‘(x+),  for  x in  (0,  »). 

n-*» 

Gn(x)  ^ ^Cx)>  f°t  n ■ 1,  2,  ....  and  every  real  number  x. 

Leirana  2.2.  Let  x and  y be  real  positive  numbers,  then 


(2.1) 

(2.2) 

(2.3) 

(2.4) 

(2.5) 

(2.6) 


K . 

/ F' (u+l 
x F (u) 


du  = JlnF  (y)  - JtnF(x). 


(2.7) 


Proof.  Since  by  (2.6)  G (x)  and  G (y)  are  positive  real  numbers,  we  obtain 
y gn(u)  n n 

that  > G '(u')du  * AnGn(y)  " ^nGn ^ * Statement  (2.7)  follows  from  (2.4),  (2.5) 
^ n 

and  the  monotone  convergence  theorem. 
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Let  £(x)  be  the  cause  of  occurence  function  given  by 

{ 1 } I (T2  > xT^)  ♦ {2}I(T2  < xTj)  ♦ {1,  2}I(T2  * xTj).  (2.8) 

Lemma  2.3.  If  £(x)  and  T(x)  are  independent  random  variables  for  every  real 
number  x,  then  for  z in  (0,  «>) 

* P(T,  uT  ) 

F(z)  « v2exnr.f  — /_ i—Cu).  (2.9) 

Proof.  By  (2.3)  and  the  independence  assumption  F’(z+)  equals  to 
z 1F(z)P(T2  > zTj) . Equation  (2.9)  is  obtained  from  (2.2)  and  (2.7). 

Theorem  2.4.  Let  £(x)  and  T(x)  be  independent  random  variables  for  every  real 
number  x.  Further  let's  assume  that  Lim  k'^ET^)1^^  is  finite.  Then 

k-H»  1 

T(x)  and  u^Ffx)^  are  identically  distributed  for  x in  (0,  ») . (2.10) 

Proof.  From  (2.9)  we  obtain  that  E Tk(x)  = ECy^FWT  ]k  for  k = 1,  2,  ...  . 
Consequently  (2.10)  follows  from  the  moments  property  of  T„.  [Breiman  (1968), 
pp.  182,  proposition  8.49]. 

Remark  2.5.  Since  x min(x'1T1,  T2)  = minO^,  xT2),  it  follows  that  if  the 
conditions  of  Theorem  (2.4)  are  satisfied,  then 

Ul^Tj  and  ^2^2  are  identically  distributed.  (2.11) 


Lemma  2.6.  Let  T^  and  T2  have  positive  continuous  densities  f^  and  f2 
respectively.  Further  let  x and  y be  positive  real  numbers,  and  a » yx”1.  If 
(2.10)  is  satisfied,  then  the  following  two  equations  hold. 
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P{T2  > y|T1  = x}  = F_1(a))a2  F'2(a)F' (a*) . (2.12) 

P{T2  * y|Tj  = x)  = f j1 (x)f2(y  F1  (a) )a2  F"2(a)F'  (a-) . (2.13) 

Proof.  Equation  (2.13)  follows  clearly  from  (2.12).  To  prove  (2.12)  it 
suffices  to  consider  the  case  F'(a+)  = F'(a-).  A proof  for  this  case  involves 
elementary  calculations  and  is  omitted. 

Remark  2.7.  Let  Z be  a positive  random  variable  independent  of  (T^,  T2).  If 
the  vector  (T.,  T2)  satisfies  (2.10),  then  so  does  (T^Z,  T2Z) . 

Theorem  2.8.  If  (2.10)  holds,  then  for  every  real  number  x the  random  varia- 
bles £(x)  and  T(x)  are  independent. 

P:~oof . Firstly  let's  assume  that  T^  and  T2  have  positive  and  continuous 
densities  f^  and  f2  respectively.  Since  PfT(x)  > t,  C(x)  = (1)1  = 

OO  00 

/ P[T2  > xu|T  = ulf  (u)du,  and  PfT(x)  > t,  £(x)  = {2}]  = / P[xu  > T2  > t|T  =u]f/u)du, 
t/x  t/x 

the  result  for  this  particular  case  follows  from  (2.12)  and  (2.13).  To  complete 
the  proof,  let  Z^,  Z2,  ...,  be  a sequence  of  random  variables  independent  of 
(Tj,  T2),  given  by  P[Zn  s t]  - [1  - e'n(t'1)ll(t  2 1),  n = 1,  2,  ...  . T^,  T^ 

have  continuous  positive  densities  and  by  Remark  (2.7)  satisfy  condition  (2.10) 
for  n = 1,  2,  ...  . Consequently  Z^Tfx)  and  C(x)  are  independent  for  every  real 
number  x,  and  positive  integer  n.  Since  Z^  converges  in  probability  to  1 as 
n -*•  » the  desired  result  follows. 


Finally  we  note  that  if  £(x)  and  T(x)  are  independent  for  every  real  number  x, 
then  by  (2.3)  the  following  tnree  equations  hold  for  z in  (0,  <») . 
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>2^1  = 2 F'  (z+)/F (z)  . 

(2.14) 

Z Trl  » 1 - 2 F'(z-)/F(z). 

(2. IS) 

] = zC(F»  Cz-)  - F1 (z+) ]/F (z) . 

(2.16) 
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3.  Main  Results. 

Let  (T  , ....  T ) be  a positive  random  vector  with  means  equal  respectively 

to  jij,  . ..,  and  let  I be  the  set  °f  all  nonempty  subsets  of  {1,  . ..,  n}. 

Further  let  T(a  , . a ) be  min  a.T. , and  let  F(a  , . ..,  a ) be  the  expected 
1 n lsisn  11  in 

value  of  T(a  , . ..,  a ).  We  define  the  cause  of  occurence  function  £(a  , a ) 


l j I (min  a.T.  < min  a T. ) . 
JeJ  ieJ  1 1 


(3.1) 


Finally  let  denote  by  statement  (3.2)  the  following  property  of  (Tj , ...,  Tn) , 

There  exists  an  integer  iQ  in  (1,  ...»  n),  such  that  T(a^ aR)  and 

T.  v . * F(a1,  . ..,  a ) are  identically  distributed  for  every  n positive  real 

*0  x0  1 n 

numbers  a. , ....  a . 

1 ’ n 

For  reference  purposes  we  note  that 


(3.2) 


Lemma  3.1.  If  (3.2)  is  satisfied,  then  the  following  three  statements  hold. 

P *T  f r = 1,  ...,  n are  identically  distributed  random  variables.  (3.3) 

For  J e I and  a.,  i e J positive  real  numbers,  min  a.T.  and  p'*(E  min  a.T.)"*T 
i . t ii  r . y l l r 

ieJ  ieJ 

are  identically  distributed  random  variables  for  r = 1,  ....  n.  I'5-  J 

For  J e J and  a.,  ...,  a positive  real  numbers,  min  a.T.  and 

1 n lsi<n  1 1 (3.5) 

(min  a.T. )F(a.,  ...,  a ) (E  min  a T ) are  identically  distributed  random  variables. 

• . ii  jl  n • . it 

ieJ  ieJ  1 1 

We  are  ready  to  extend  Theorems  2.4  and  2.7. 


Theorem  3.2.  Let  ?(a^,  ....  a^)  and  T(a^,  ...,  a^)  be  independent  random  variables 


for  every  n real  numbers  a,  , . . . , a . Further  let  assume  that  Lim  k T ) 

I n . r' 

k-*® 

is  finite  for  some  positive  integer  r in  {1,  ....  n).  Then  (3.2)  holds. 


k.l/k 
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_ i i,  i/i. 

Proof . Let  Lim  k (E  T.  ) be  finite,  and  let  U and  U be  respectively 
k-*=°  0 1 2 

equal  to  min  a.T.  and  a.  T . Since  (U. , U_)  satisfies  the  conditions  of 
i*iQ  X0  X0 

Theorem  (2.4),  statement  (3.2)  follows. 

Theorem  3.3.  If  (3.2)  holds,  then  for  every  n real  numbers  a a . the 

I n 

random  variables  C(a^»  • a^)  and  T(a^,  ...,  a^)  are  independent. 

Proof.  Let  Jel  and  be  equal  respectively  to  min  a.T.  and  min  a.T. . 

ic.T  11  i^J  1 1 

Since  by  (3.5)  condition  (2.10)  is  satisfied,  the  result  is  obtained  by 
Theorem  2.7. 

Remark  3.4.  Let  J be  in  I,  further  let  , U7  be  respectively  equal  to 

min  a T and  to  min  a T. . Finally  let  H(x)  be  equal  to  E nin(xU. , U_) . If 
ieJ  11  i^J  1 1 12 

(Tj,  ....  T ) satisfies  (3.2),  then  by  (2.3) 

Pts(aj,  ••.,  a^)  = J)  = H' (x+) lx=j  * F (a^ , ...,  a^).  (3.6) 

Remark  3.5.  Let  Z be  a positive  random  variable,  independent  of  (T^ , ...,  T ). 

If  (Tj,  ....  Tn)  satisfies  (3.2),  tlien  so  does  (T1Z,  ....  T Z) . 

A nonnegative  random  variable  T has  a Weibull  distribution  with  parameters 
p and  n if  for  t in  (0,  «) 

P(T  > t}  = exp[-pta],  p,  a > 0.  (3.7) 

Let  a be  a positive  real  number  and  g a positive  real  function  defined  on  the 
nth  Euclidian  space.  !Ve  say  that  a nonnegative  random  vector  (T, , ....  T ) has 
Weibull  minima,  if  for  every  n positive  real  numbers  ....  a^  and  t in  (0,  «°) 

PtTfaj,  •••»  an)  > t)  = exp[-g(a1,  ....  a^Jt0].  (3.8) 


9 


Clearly  every  subset  of  n random  variables  with  Weibull  minima  has  Weibull 
minima.  If  T is  a Weibull  random  variable  with  parameters  u and  a then 

oo 

_ _k  -k/a  r -z  k/a,  ,, 

E T = vi  J e z dz.  (3.9) 

0 

Consequently  if  T is  Weibull  with  a £ 1 then 

Lim  k-1(E  Tk)1/k  < >*>  (3.10) 

k+<» 

and  the  following  corollary  holds: 

Corollary  3.6.  (i)  If  (1^,  ...,  T^)  has  Weibull  minima  then  5(aj,  ...,  afl)  and 

T(a^,  . ..,  a^)  are  independent  random  variables  for  every  n positive  real 
numbers  a^,  ....  an>  (ii)  If  at  least  one  of  the  T^'s  is  Weibull  with  a 2 1, 
and  C(a^,  . ..,  a^)  is  independent  of  T(a^,  ....  a^)  for  every  n real  numbers 
al’  *•*»  an’  t*13n  • •.,  Tn)  has  Weibull  minima. 

Proof.  Follows  from  Theorems  3.2,  3.3,  and  from  (3.10). 

Essary  and  Marshall  (1974)  defined  a class  of  positive  n-dimensional  random 

vectors  that  have  the  following  property. 

For  J e J and  a. , i e J positive  real  numbers,  min  a.T.  has  an  exponential 

ieJ  1 1 (3.11) 

distribution. 

Corollary  3.6  provides  in  particular  the  following  characterixation  of  that 
class. 

Corollary  3.7.  Let  (Tj,  ....  T^)  be  a positive  random  vector  with  at  least  one 
exponential  component.  (Tj,  ....  Tn)  satisfies  (3.11)  iff  £(aj,  ...,  an)  and 
T(a^,  . ..,  a^)  are  independent  random  variables  for  every  n real  numbers  a^,  ...,  an. 
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Let  R(t)  be  a positive,  strictly  increasing  function,  that  is  differentiable 

and  converges  to  ® as  t ■*■  «.  Further  let  h be  a positive  real  function  on  the 

nth  Euclidian  space.  We  say  that  (Tj , T^)  has  proportional  minima,  if 

for  every  n positive  real  numbers  a,,  ....  a and  t in  (0,  =°) 

i n J 

P{T(ai,  ....  an)  > t}  = expr-hfa^  ....  an)R(t)].  (3.12) 

Finally  we  show  that  the  only  family  of  random  vectors  with  proportional  minima 
is  the  Weibull  one. 

Theorem  3.10.  Let  (T^,  ...,  T^)  be  a positive  random  vector.  (T^,  ...,  T ) (n  ^ 2) 
has  proportional  minima  iff  it  has  Weibull  minima. 


Proof.  Let  6(a)  be  h(a,  ...,  a),  then  R(ta)0(a)  = R(t)0(l).  Consequently 
t dR  _ a dh  , . . . c 

R dt  = h da’  lience  = At  for  some  positive  real  numbers  A and  c. 
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4.  Examples. 

Example  4.1.  Let  A ,,  J e J be  nonnegative  real  numbers,  that  add  up  over  the 
sets  in  J to  a positive  number.  We  define  the  multivariate  distribution  of 
(Tj,  ....  Tn)  for  tj,  ....  tn  in  (0,  »)  by 

P{T.  > t.,  i = 1,  . ..,  n)  = expC-  £ A max  t“],  a > 0.  (4.1) 

1 JeIJieJ  1 

The  distribution  in  (4.1)  is  an  extension  of  the  bivariate  Marshall-Olkin 

(1967)  exponential  distribution.  Clearly  this  family  of  random  vectors  has 

Weibull  minima  with  g(a.,  ....  a ) = T A. max  aTa.  In  particular  we  obtain 

J el  leJ 

for  Aj  * 0,  whenever  J is  not  a singlton  set  independent  Weibull  random  variables. 

If  Aj  = 0 for  J * (1,  ...,  n}  the  multivariate  distribution  given  by  (4.1) 
reduces  to  correlated  Weibull  random  variables. 

Example  4.2.  Let  Aj,  A2  be  nonnegative  real  numbers.  Further  let  Yj,  Y2  ••• 
be  a sequence  of  positive  numbers  that  add  up  to  a finite  real  number.  Finally 
let  bj,  b2  ...  be  a sequence  of  positive  real  numbers.  We  define  the  bivariate 
distribution  (T^,  1^)  for  tj,  t2  and  0 in  (0,  <*>)  by 

CO 

P(T.  > t.,  i * 1,  2}  = expC-Ajt ° - A2t“  - [ A max(t“  tV“)].  (4.2) 

n=l 

This  distribution  has  Weibull  minima  with  g(a1>  a ) equal  to 
00 

A.a"a  ♦ A_a~a  ♦ £ A (a"aa _°b’a) . One  can  extend  with  no  difficulty  (4.2)  to 

l l n«i  n 1 ^ n 

the  multivariate  case,  although  the  explicit  expression  becomes  some  what 
cumbersome . 


Using  Remark  (3.5)  we  can  generate  the  following  two  examples. 
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Example  4.3.  Let  (T^ , T^)  be  given  by  (4.1)  and  let  Z be  a Weibull  random 

variable  with  parameters  a and  a independent  of  (T^ , . . . , T ) . The  random 
vector  (Z  *T^,  ....  Z *T^)  with  a survival  probability  at  tj,  ....  t^, 

(min  tj  > 0) , given  by 

yatn  + l A max  ta]  1 (4.3) 

■T  J iej 

satisfies  (3.2). 

Example  4.4.  Let  (T^,  . ..,  T ) be  given  by  (4.1)  and  let  Z be  a positive  random 

variable  independent  of  (T  , . . . , T ) with  a density  equal  to 

11 

e wt  t®  ^ap  a T^e  random  vector  (Z_1T^,  . ..,  Z-1^)  with  a 

survival  probability  at  tj,  ....  t , (min  t.  > 0)  equal  to 

M,i 

Tu(u  + £ X T max  ta)  a (4*4) 

J ieJ  “ 

satisfies  (3.2). 

Finally  let  Pj,  . . . , be  positive  real  numbers,  let  (3  be  in  the  set 

(0,  1]  and  let  a be  in  [1,  <*>) . The  survival  probability  given  by 

expC-  l p.t“]B,  (4.5) 

i=l  1 1 

n 

clearly  has  Weibull  minima  with  g(a  , ...,  a ) equal  to  ( J p.a7°)  . 

i=l  1 1 

Essary  and  Marshall  (1974)  presented  two  bivariate  distributions  that  have 
exponential  marginals  and  the  minima  of  the  two  components  is  exponential,  .o’-'ever 
their  joint  distributions  are  not  bivariate  Marshall-Olkin  (1967)  exponential. 

The  absolute  continuous  distribution  given  in  (4.5)  for  <?.  = 1/a  is  a si  iple 
r.-dimensional  example  for  such  a situation. 
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